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We use the pubhcly available code CMBFAST as modified 
by L. Pogosian and T. Vachaspati to simulate the effects of 
wiggly cosmic strings on the Cosmic Microwave Background 
(CMB). Using the modified CMBFAST code which takes into 
account vector modes and models wiggly cosmic strings by 
the one scale model, we go beyond the angular power spec- 
trum to construct CMB temperature maps with resolution of 
a few degrees. The statistics of these maps are then stud- 
ied using conventional and recently proposed statistical tests 
optimized for the detection of hidden temperature disconti- 
nuities induced by the Gott-Kaiser-Stebbins effect. We show 
however, that these realistic maps can not be distinguished 
in a statistically significant way from purely Gaussian maps 
with identical power spectrum. 

PACS:11.27,98.80C,02.50,07.57H 



I. INTRODUCTION 

A promising method of identifying the origin of the 
large-scale structure of the universe, is the statistical 
analysis of the CMB fluctuations [1-5]. What one hopes 
to get out of this is the possible distinction between the 
two major classes of physical theories of structure for- 
mation, namely theories of inflation [pi and topological 
defects [0 . An important differentiating feature of these 
two classes of theories is the statistics of primordial fluc- 
tuations. Temperature anisotropics -^ of the CMB, if 
originating from quantum fluctuations in an inflationary 
context, should obey Gaussian statistics. However, the 
confidence by which the CMB anisotropics were shown 
to follow that statistic is reduced M , taking in mind the 
work done by §,|l§ using COBE-DMR maps. This could 
favor topological defect theories. 

Cosmic strings seeding structure formation [11-15] con- 
sist a promising case of these theories. However, despite 
the significant progress made recently [12,13,16-18] in un- 
derstanding the effects of cosmic strings on the CMB, 
the uncertainty of the derived predictions remains sig- 
nificant ll^,p^. The main source of this uncertainty is 
that there is no simple way to characterize the network of 
strings. In order to bypass this difficulty various studies 
have attempted to model the string network by incor- 



porating simplifying assumptions. Such simplifications 
have made possible the realization of cosmological string 
network simulations [21-23] that attempt to model evo- 
lution in a very wide range of scales. They have also 
made possible the construction of semi-analytical models 
that attempt to capture features missed by the inherent 
limitations of full scale numerical simulations. For exam- 
ple, early work [p4,pq| based on the latter approach had 
revealed important features of the cosmic string induced 
CMB fluctuations that were missed by detailed numeri- 
cal simulations. In particular it was shown [ p5[ , using a 
simple model of string network evolution, that the wig- 
gles (small scale structure) of long cosmic strings tend 
to amplify the height of the Doppler peak of the CMB 
angular spectrum of string induced fluctuations. 

The main reason for this amplification may be at- 
tributed to the two main factors that determine the 
string induced fluctuations. These are, first the Inte- 
grated Sachs-Wolfe effect (long strings present between 
last scattering and the present time interacting gravita- 
tionally directly with the photons) for scales larger than 2 
degrees and second, for scales about 1-2 degrees, interac- 
tions of the photons with the string perturbed plasma 
performing acoustic oscillations on the last scattering 
surface (LSS) excited by long strings. It may be shown 
[25-27] that the latter effect leads to temperature fluctua- 
tions =p^ that are amplified compared to the fluctuations 
induced by the former effect. The amplification factor 
A > 1 is due to the fact that the wiggles do not interact 
gravitationally directly with the photons, while they do 
interact with the plasma [£5|,^ on the LSS leading in- 
directly to additional fluctuations on the CMB through 
the last scattering of photons on the fast moving plasma 
electrons. Therefore this amplification can only affect 
the last scattering horizon scale, i.e. a scale of about 1-2 
degrees. The factor A can be expressed as 25 27[| 



A = (l + 
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(1.1) 



where T is the tension of the wiggly long string which 
decreases with the wiggliness and is estimated by sim- 
ulations to be T ~ 0.7/i where /i is the mass per unit 
length of the string. Also, Vs is the string velocity and 
7s the corresponding Lorentz factor. 



This result of Doppler peak amplification due to cosmic 
string wiggles has recently been verified by the numeri- 
cal modeling developed by Pogosian and Vachaspati [g6| . 
These authors used the line of sight integration approach 
of the publicly available code CMBFAST (2| combined 
with the one-scale model [^ 30 1 for wiggly cosmic string 
evolution to derive the angular CMB power spectrum and 
the matter power spectrum. One of their main results 
was that small scale structure (wiggles) of long strings 
tends to increase the height of the Doppler peak of the 
CMB spectrum thus confirming the expectations of Ref. 
p5[ and improving the agreement of the predicted angu- 
lar spectrum with observations. 

In the present work we use the approach of Pogosian 
and Vachaspati and modify their publicly available code 
based on CMBFAST to go beyond the CMB spectrum 
and construct realistic maps of CMB fluctuations pro- 
duced by wiggly strings. Our goal is to apply statisti- 
cal tests on several realizations of these maps and at- 
tempt to identify their non-Gaussian features. Some of 
the tests we applied |3l) (which were specially designed 
for this purpose) were able to identify non-Gaussian fea- 
tures on simple CMB maps Q , constructed by the super- 
position of a small temperature discontinuity (produced 
by a single long string via the Gott-Kaiser-Stebbins effect 
|p2| , p3| ) on a Gaussian background. In the present study 
however, we show that these same tests are not able to 
identify the string induced non-Gaussianity of more re- 
alistic maps on the scales considered. This is due to the 
fact that the superposition of several line-like tempera- 
ture discontinuities tends to wash out the effects of each 
individual discontinuity and reduces the effectiveness of 
the proposed tests. Therefore, the detection of the small 
non-Gaussianity induced by cosmic strings on a few de- 
grees scale remains a challenging open issue. 



II. CONSTRUCTION OF THE MAPS 

The construction of the CMB temperature anisotropy 
maps on which the statistical tests were performed, was 
made possible by two major contributions: the wiggly 
string one-scale model as the cosmic string network 
model and the use of the line of sight integration ap- 
proach. 

The cosmic string network model that we used is the 
one found in detail in p6[ . Briefly, the network is approx- 
imated [12-14] by a collection of uncorrelated straight 
string segments which are assumed to be produced at 
some early epoch, with positions drawn randomly from a 
uniform distribution in space, moving with random un- 
correlated velocities, their directions drawn from a uni- 
form distribution on a two sphere. This model has the 
merit of being relatively simple and amenable to modi- 
fications that seem to be indicated by direct simulations 

Small-scale structure, wiggles, were incorporated in or- 



der to make the model more realistic, since lattice sim- 
ulations of string formation and evolution |g^ suggest 
that strings are not straight. However, a distant observer 
would not be able to discern this small-scale structure, 
but see instead a smooth string with effective mass per 
unit length U and tension T. The wiggly string is heav- 
ier and slower than the ordinary Nambu-Goto straight 
string (whose equation of state is: hq = T, where /io is 
the mass per unit length). The equation of state of a 
wiggly string (averaged over the small-scale structure) is 
p5|,p6|: 



Uf = Mo, 
U = a^o, T = fio/a 



(2.1) 



where a is the "wiggliness" parameter, estimated l23] 
in the radiation and matter eras to be a^ — 1.9 and 
a,„ ~ 1.5 respectively. The expected evolution of the 



(2.2) 



wigghness parameter is fitted by |23|J2 



a(r) = 1 + 



(ofr — l)a 



where a is the scale factor. 

The parameters of the segments: length, velocity 
and wiggliness are modeled using the one-scale model 
p9| , p0[ . The evolution of the two competing processes of 
string stretching on the one hand (due to cosmic expan- 
sion) , and the chopping off of loops and their subsequent 
decay on the other (due to long string reconnection) , are 
described by [B7|: 



dl a 2 1 ~ 

—- = —Iv + —cv 
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(2.3) 
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where I is the comoving correlation length, v is the rms 
string velocity, c is the loop chopping efficiency and k 
is the effective curvature of the strings [Of. At every 
subsequent epoch, a certain fraction of the number of 
segments decays in a way that preserves scaling. 

The Fourier transform of the energy-momentum tensor 
for an individual string segment m is: 
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while 0JJJ can be found from the relation: V^O^j, = 0. 
(A^(cr, r) are the coordinates of the segments: X'^ = 
T, X = xo+crX'+wrX, where xq is the random location 

of the center of mass, X' and X are unit vectors along and 
perpendicular to the string which are randomly oriented, 

also satisfying X' • X = and a is the coordinate along 
the string). 



A consolidation of all string segments that decay at 
the same epoch into one with the appropriate statistical 
weight is used, as was suggested in Q, so that the num- 
ber of segments can be dealt with computationally. The 
total energy-momentum tensor of the network is the sum 
over the contributions of the consolidated string segments 



where the term S{k, r) contains the contributions from 

all the sources, as well as Q^^, from the string network 
p8| . Using the relation: 



Jkq(T-Ta) 



J2i-^y{'2l + m[kiro-T)]Pi{q) (2.13) 
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e^.(k,T) = J2 0;'l(k,r)r°/^(T,T™), 



(2.7) 



where Nq is the initial number of segments, and T"^^ 
is a smooth function that turns off the ttt.*'' segment by 
time Tm ||l^ , |l3|] . 9my(k, t) will be incorporated into the 
sources: S{k,T) (cf eqs.( |2J^ , |2l^ )). 

Our main goal is to compute the CMB anisotropy 
seeded by cosmic strings: 
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(n) 
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(x = 0, n,To) =^aimyi™(n) (2.. 



(x is the position of the observer and n is the line of 
sight direction). The codeQof |§ (based on CMBFAST 
m which we modified and used, uses the line of sight in- 
tegration approach [g8[ to compute, among other things, 
the angular power spectrum Cf. 



1 ' 
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(2.9) 



originating from temperature perturbations seeded by 
cosmic strings by incorporating the sources into 0^^ in 
the way described below. 

The Fourier transform A(k, n, ro) of ^T^(n) depends 
only on the angle q between the two vectors (q = k-n) and 
can be expanded in multipole moments A/(k, tq) ||38|| : 

CXj 

A(k,n,To) = ^(-z)'(2; + l)A,(k,To)P/(k-n) (2.10) 

1=0 

further, A/(k, tq) is decomposed into Q: 



Az(k,To) = ^(k)A,(fc,To) 



(2.11) 



However, according to the line of sight integration ap- 
proach |2a,p3]: 



A{k,q,To)= e''"i(^-^°^S{k,T)dT (2.12) 

^0 



we get 



AKfc,To)=/ Sik,T)jl[k{To-T)]dT (2.14) 

Jo 

In order to get the aim coefficients needed to construct 
the maps, we have to assume arbitrary direction for k, 
or equivalently use the more general relation n. 

oo / 

e^kr = 4nJ2Yl ^'jiikr)Yi:,ik)Yim{r) (2.15) 

/=0 m=-l 



in ( 2.10| ,2.12). Otherwise we would be able t o com pute 
only the aio terms as we shall see next. Then (2.14) still 
holds. The aimS read PSJ]: 



ai,j 
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dnY^*min) — {n)=A7^{-^y / d^kYC^MM^To) 



Ani^tY /d3kr;:jk)e(k)A,(fc,To) 



(2.16) 



where A/(fc,ro) is given in (2.14) (had we not assumed 
the arbitrariness of k, we would only get: aim — 
{-iy^M'2l + l)JdM{k)Ai{k,To)Smo^aio). 

Since the linearized Einstein equations may be decom- 
posed into scalar, vector and tensor parts (S, V, T) Q, 
we have three distinct contributions to the sources (e.g. 
S: density perturbations, V: cosmic strings, T: gravi- 
tational radiation, [OJlql). So, (A, A;) are actually: 
(A^, Af ), where I = S,V, T. Then: 
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(2.17) 



while Ci = Ci + CY + Cf . Following the normalization 
used in the code, the a/mS are defined as: 



- / d^ky,:„(k)Af(k,To) 



''Im 



''Im 



1(1 + 1) d^kYC,,{k)AY{k,To) (2.18) 



U {i + 2y. 

n (1-2)1 



d3ky4(k)Af(k,To) 



where: A[(k, tq) = ^(k)A^(fc, tq). In order to implement 
this scheme, we had to extend the code of [E6| to provide 



^http://theory4. phys.cwru.edu/~levon 
^http://www. sns.ias.edu/~matiasz/CMBFAST/cmbfast 
.html 



L.D. Landau, in Quantum Mechanics, Chap. V, § 34, p. 



113. 



for the amplitude of the initial perturbations £,(k), which 
must be such that IsslMEll : 



<r(k)C(k')>=^«(fc)^(k-k') 



(2.19) 



where P^(k) is their power spectrum. 

As in ||26|] , we chose the initial power spectrum to be 
that of white-noise, i.e. P^{k) = constant. Consequently, 
^(k) was chosen to be proportional to a uniformly ran- 
dom number: ^k e[~lil]- Since: < ^kCk' >= l^kk', each 
'A.m. ^^^ multiplied by a factor of vS when the ensemble 
average < aj*^aj„^ > was taken in order to compute and 
verify the C/s Q 

Applying: J d^k' < ^*(k) C(k') >= F^(fc) to our dis- 
crete case we get: Y.k ^^k < Ck^k' >= 5 Lk ^^kf^kk' = 
|AVk. This leads to the definition of ^(k) as: 



e(k) 



a 



VAK 



^k '■ random e[— 1, 1] 



AVk = k^dk dn(k) 



(2.20) 

Having chosen this particular realization for ^(k), we 
can proceed to the const ruction of the temperature 
anisotropy maps using eqs.(2.8, 2.18 ). 



III. RESULTS 

We chose to construct maps with / ~ 40, so that the 
pixel size is comparable to that of the COBE maps, al- 
though our method enables us to proceed even up to 
much higher Is (e.g. I ~ 1500), with greater cost in com- 
putational time though. The angular power spectrum of 
such a map is depicted in Fig.n^ in agreement with Ref. 




FIG. 1. Angular power spectrum of a Z = 1500 map. 



*The C/s are: Cf = ^ J Pdk\Af{k,ro)\'' , CY = ^l{l 



A typical example of such a map is depicted in Fig||, 
where a 30° x 30° patch of a standardized temperature 
anisotropy map is shown. 
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FIG. 2. A 30° X 30° patch of one of the constructed stan- 
dardized temperature anisotropy maps. (Lighter colors cor- 
respond to higher values of ^-). 

We have applied statistical tests on the above maps 
in an attempt to identify possible non-Gaussian features. 
Such statistical tests [Q have recently been proposed for 
the detection of the particular type of non-Gaussianity 
induced by coherent discontinuities (edges) hidden in 
CMB maps. The main advantage of these statistical 
tests is that they focus on the large-scale coherent prop- 
erties of CMB maps and are therefore effective even in 
cases of low resolution maps provided that the area cov- 
ered is large. The statistics calculated for maps of CMB 
anisotropics seeded by cosmic strings which we have con- 
structed are the skewness, kurtosis and the recently de- 
vised MSD, SMD @ which are optimized for the de- 
tection of coherent discontinuities in 1-D and 2-D pixel 
maps. It should be pointed out that the skewness and 
kurtosis statistics have proven pl| to be notably unsuit- 
able to detect the subtle non-Gaussianity in CMB maps. 
However they have been used as a benchmark [ pl[ to 
test the effectiveness of more complex statistics like the 
SMD and the MSD specially designed for the detection 
of specific types of non-gaussian signatures. 

In addition to the statistics discussed here, there are a 
number of other more complex statistics p2] (minkowski 
functionals and bispectrum and higher order cumulant 
estimates in harmonic space) that may be applied to 
CMB maps to check for non-Gaussianity or topological 
signatures of defects. Those statistics which are partic- 
ularly effective on high resolution CMB maps have been 
discussed extensively in the literature and their applica- 
tion on string maps produced by CMBFAST can be an 
interesting extension of the present work. Here we focus 



on the MSD and SMD statistics which are specially de- 
signed P to detect the specific signature of large scale 
coherent discontinuities on CMB maps of low resolution. 
Assuming a 30° x 30° (in pixels) standardized tem- 
perature map Tij{i,i — 1,...,30), the skewness s and 
kurtosis k are defined as : 






(3.1) 

(3.2) 



«j 



These are conventional statistics and their values for 
Gaussian maps with uncorrelated pixels are s — Q^k — i. 
To define the statistics MSD and SMD we consider a 
partition of the CMB maps in two parts separated by a 
random curve. In this study we have considered straight 
lines as well as right angles. Let k denote the set of pa- 
rameters that define the partition line and let T^ and T; 
be the mean temperatures of the two parts of the map 
(the indices "u" and "1" stand for "upper" and "lower" 
parts). The statistical variable Yj. is defined as : 



n = \Tu 



T,\ 



(3.3) 



The statistics MSD and SMD are defined as: 



N 



SMD = ia&yi{Yj?), 



(3.5) 



where N is the total number of partitions. Both MSD 
and SMD approach asymptotic values for large- A^. Con- 
sider now a Gaussian pattern of temperature fluctuations 
with a small coherent temperature discontinuity defined 
by a partition /cq superposed on the map with coherence 
scale comparable to the size of the pattern. In ||31|] it was 
shown that the presence of such a discontinuity can be de- 
tected much more efficiently by the statistics MSD and 
SMD than by the skewness and kurtosis. A physically 
motivated mechanism which can lead to the production 
of a coherent discontinuity on CMB maps is the presence 
of a moving long cosmic string in our horizon via the 
fore- mentioned Gott-Kaiser-Stebbins effect p2| , p3[ . 

The four statistical tests (MSD, SMD, skewness, kurto- 
sis) have been performed on patches of the maps (Fig.||) 
produced by the simulation in a way similar to the one 
in§: 

• First the values of the 4 statistical variables were 
computed for a specific 30° x 30° patch obtained 
using the simulation code. 

• A large number (1200) of Gaussianized maps were 
constructed from this patch by randomizing the 
phases in Fourier space and using a Gaussian spec- 
trum with mean and variance equal to the mea- 
sured spectrum of the map. 



• The 4 statistics were calculated for each of these 
Gaussianized maps. 

• Using these results, the probability distribution for 
each statistic was constructed. 

• Finally, the probability for obtaining the already 
calculated (in step 1) values of the 4 statistics was 
found. 

The results are shown in Figs. || |j. Clearly, the ap- 
plied statistics can not reveal non-Gaussianity in a statis- 
tically significant amount in the maps produced by the 
simulation. The values of the statistics considered for 
the string induced maps (black dots) are within la in 
the probability distributions produced by the gaussian- 
ized maps. 
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FIG. 3. Probability distribution of the MSD based on 1200 
Gaussianized maps. The probability for obtaining the MSD 
value of the string-induced map (MSD=7.378, represented by 
the dot) is 0.35. 
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FIG. 4. Probability distribution of the kurtosis based on 
1200 Gaussianized maps. The probability for obtaining the 
kurtosis value of the string-induced map (kurtosis=3.08, rep- 
resented by the dot) is 0.17. 

We have also extended the above results using right 
angles as partition lines instead of straight lines. How- 
ever, the sensitivity of the statistics considered was not 
increased, rendering results similar to Figs, y, |4|. 



IV. CONCLUDING REMARKS 

String-induced maps on large angular scales are hard 
to distinguish from maps with Gaussian fluctuations, 
even using specially designed tests. We have suggested 
a technique for constructing string-induced temperature 
anisotropy maps which can be exploited in many ways. 
The maps constructed the way described in section 2 are 
proved not to be dominated by late long strings. On 
the contrary, small-scale structure features come into 
play which, in a way, destructively interfere and ruin the 
large-scale coherence which the MSD and SMD tests are 
proposed to optimally detect. Thus, we offered a con- 
firmation to the fact that string-induced temperature 
anisotropy maps are practically Gaussian on large scales. 



V. ACKNOWLEDGEMENTS 

We would like to thank Lcvon Pogosian and Tanmay 
Vachaspati for their assistance in acquiring and better 
understanding the code they developed by modifying 
CMBFAST. 



[1] J. Gott ei al, Ap. J 352, 1 (1990). 

[2] R. Moessner, L. Perivolaropoulos, and R. Brandenberger, 

Ap. J. 425, 365 (1994). 
[3] J. Magueijo, Phys. Rev. D 52, 4361 (1995). 
[4] P. G. Ferreira and J. Magueijo, Phys. Rev. D 56, 4578 

(1997). 
[5] A. Lewin, A. Albrecht, and J. Magueijo, Mon. Not. R. 

Astron. Soc. 302, 131 (1999). 
[6] D. Lyth, Phys. Rept. 314, 1 (1999). 
[7] A. Vilenkin and E. Shellard, in Cosmic Strings and other 

Topological Defects, edited by P. Landshoff (Cambridge 

University Press, Cambridge CB2 IRP, 1994). 
[8] A. Kogut et al, Ap. J. Lett. 464, L29 (1996). 
[9] L. Perivolaropoulos and N. Simatos, Mod. Phys. Lett. A 

13, 2945 (1998). 
[10] P. G. Ferreira, J. Magueijo, and K. Gorski, Ap. J. 503, 

LI (1998). 
[11] U. Pen, U. Seljak, and N. Turok, Phys. Rev. Lett. 79, 

1611 (1997). 
[12] A. Albrecht, R. Battye, and J. Robinson, Phys. Rev. Lett. 

79, 4736 (1997). 

[13] A. Albrecht, R. Battye, and J. Robinson, Phys. Rev. D 

59, 023508 (1999). 
[14] R. Battye, J. Robinson, and A. Albrecht, Phys. Rev. Lett. 

80, 4847 (1998). 

[15] C. Contaldi, M. Hindmarsh, and J. Magueijo, Phys. Rev. 

Lett. 82, 679 (1999). 
[16] N. Turok, U. Pen, and U. Seljak, Phys. Rev. D 58, 023506 

(1998). 



[17 

[is; 

[19 

[20 

[21 
[22' 
[23^ 

[24' 
[25^ 
[26; 

[27 

[28' 
[29^ 
[30 
[31 
[32' 
[33 
[34; 

[35' 
[36 

[37 

[38' 
[39 

[4o; 

[41 

[42; 



U. Seljak, U. Pen, and N. Turok, Phys. Rev. Lett. 79, 
1611 (1997). 

U. Seljak, U. Pen, and N. Turok, Phys. Rev. Lett. 79, 
1615 (1997). 



A. Riazuelo, N. Deruelle, and P. Peter, |astro-ph/9910290 
(1999). 

R. Durer and M. Sakelariadou, Phys. Rev. D 56, 4480 
(1997). 

B. Allen et al, Phys. Rev. Lett. 79, 2624 (1997). 
B. Allen et al, Phys. Rev. Lett. 77, 3061 (1996). 

D. Bennett and F. Bouchet, Phys. Rev. D 41, 2408 

(1990). 

L. Perivolaropoulos, Phys. Lett. B 298, 305 (1993). 

L. Perivolaropoulos, Ap. J. 451, 429 (1995). 

L. Pogosian and T. Vachaspati, Phys. Rev. D 60, 083504 

(1999). 

T. Vachaspati and A. Vilenkin, Phys. Rev. Lett. 67, 1057 

(1991). 

U. Seljak and M. Zaidarriaga, Ap. J. 469, 437 (1996). 

T. Kibble, Nucl. Phys. B 252, 227 (1985). 

T. Kibble, Nucl. Phys. B 261, 750 (1986). 

L. Perivolaropoulos, Phys. Rev. D 58, 103507 (1998). 

J. Gott, Ap. J. 288, 422 (1985). 

N. Kaiser and A. Stebbins, JVature 310, 391 (1984). 

D. Bennett, F. Bouchet, and A. Stebbins, JVa- 

ture(London) 335, 410 (1988). 

B. Carter, Phys. Rev. D 41, 3869 (1990). 
A. Vilenkin, Phys. Rev. D 41, 3038 (1990). 

C. Martins and E. Shellard, Phys. Rev. D 54, 2535 
(1996). 

C. P. Ma and E. Bertschinger, Ap. J. 455, 7 (1995). 

U. Seljak, Ap. J. 482, 6 (1997). 

R. Crittenden and N. Turok, Phys. Rev. Lett. 75, 2642 

(1995). 

M. Zaidarriaga and U. Seljak, Phys. Rev. D 55, 1830 

(1997). 

J.Magueijo, Phys. Lett. B 342 32-39 (1995), ERRATUM- 

ibid. B352 499 (1995) 

P.Ferreira and J.Magueijo, Phys. Rev. D 55 3358 (1997) 

P. Ferreira, J. Magueijo, and J. Silk, Phys. Rev. D56 , 

4592-4603 (1997). 



